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Within  the  framework  of  boundary  conformal  field  theory,  we  evaluate  the  conductance  of  stable  fixed  points 
of  junctions  of  two  and  three  quantum  wires  with  different  Luttinger  parameters.  For  two  wires,  the  physical 
properties  are  governed  by  a  single  effective  Luttinger  parameter  for  each  of  the  charge  and  spin  sectors.  We 
present  numerical  density-matrix-renormalization-group  calculations  of  the  conductance  of  a  junction  of  two 
chains  of  interacting  spinless  fermions  with  different  interaction  strengths,  obtained  using  a  recently  developed 
method  [Phys.  Rev.  Lett.  105,  226803  (2010)].  The  numerical  results  show  very  good  agreement  with  the 
analytical  predictions.  For  three  spinless  wires  (i.e.,  a  Y  junction)  we  analytically  determine  the  full  phase 
diagram  and  compute  all  fixed-point  conductances  as  a  function  of  the  three  Luttinger  parameters. 
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I.  INTRODUCTION 

Conducting  quantum  wires,  at  low  energies,  generically 
form  a  Tomonaga-Luttinger  liquid  (TLL),  characterized  by  a 
Luttinger  parameter  g ,  which  encodes  the  effects  of  electron- 
electron  interactions.1-4  Due  to  the  prominent  role  of  interac¬ 
tions  in  low  dimensions,  the  nature  of  these  one-dimensional 
electronic  systems  is  dramatically  different  from  their  higher¬ 
dimensional  counterparts  described  by  Landau’s  Fermi-liquid 
theory.5  The  TLL  state  of  matter  in  one-dimensional  quantum 
wires  has  been  realized  in  numerous  experiments  over  the  last 
few  years.6-13 

Transport  properties  of  such  quantum  wires  are  of  consider¬ 
able  interest:  From  a  fundamental  point  of  view,  a  large  number 
of  interesting  phenomena  have  been  predicted  and  observed. 
For  instance,  at  low  temperature  and  low  bias  voltage,  a  TLL 
with  repulsive  interactions  (g  <  1)  is  totally  disconnected  in 
the  presence  of  an  impurity,  while  one  with  attractive  interac¬ 
tions  (g  >  1)  conducts  as  in  the  absence  of  the  impurity.14-19 
From  a  practical  viewpoint,  junctions  of  TLL  wires  serve  as 
important  building  blocks  of  quantum  circuits,20,21  and  are 
thus  of  technological  significance.  Junctions  of  three  quantum 
wires,  known  as  Y  junctions,  also  have  highly  nontrivial 
transport  properties.21,22  Due  to  their  rich  transport  behavior, 
junctions  of  quantum  wires  and  their  networks  have  thus 
attracted  much  attention.23-37 

Most  of  the  previous  works  on  the  transport  properties  of 
junctions  of  TLL  wires  focus  on  wires  with  the  same  Luttinger 
parameter.  However,  experimentally,  there  is  no  reason  for  all 
the  TLLs  emanating  from  a  junction  to  be  identical.  Moreover, 
a  single  TLL  can  have  inhomogeneities;  for  example,  a 
contact  between  an  interacting  TLL  and  a  Fermi-liquid  lead, 
a  key  ingredient  of  most  transport  measurements,  is  often 
studied  as  an  inhomogeneous  TLL  wire  smoothly  interpolating 
between  interacting  (TLL)  and  noninteracting  (Fermi-liquid) 
regions  or  as  a  two- wire  junction  with  the  Luttinger  parameter 
abruptly  changing  at  the  junction.38-56  A  junction  of  three 
quantum  wires  with  different  Luttinger  parameters  has  been 


studied  in  the  weak  coupling  regime.21,57-59  The  experimental 
importance  of  junctions  of  TLL  wires  with  generally  unequal 
Luttinger  parameters  motivates  an  in-depth  study  of  their 
properties,  which  is  the  main  objective  of  the  present  paper. 

Here,  we  focus  on  junctions  of  two  and  three  nonchiral 
Luttinger  liquids  schematically  depicted  in  Fig.  1.  For  two 
wires,  it  is  known  that  the  transport  properties  of  the  junction 
are  fully  controlled  by  one  effective  Luttinger  parameter 
ge  =  2/(gfl  +  gfl)  as  found  in  Ref.  39.  In  the  context  of 
fractional  Hall  edge  states,60  similar  results  have  been  found 
for  tunneling  between  two  chiral-TLL  edge  states  43  For  two 
nonchiral  wires,  one  can  reach  the  same  conclusion  through  an 
almost  identical  argument.  In  this  work,  however,  we  obtain 
this  result  within  the  framework  of  boundary  conformal  field 
theory  (BCFT),  using  the  delayed  evaluation  of  boundary 
conditions  (DEBC)  method,61-64  which,  as  we  will  see,  has 
the  advantage  that  it  can  be  readily  generalized  to  junctions  of 
more  quantum  wires. 

Such  generalization  to  a  junction  of  three  nonidentical 
quantum  wires  is  a  key  result  of  this  paper.  We  find  the  stability 
regions  of  the  previously  identified  (in  Ref.  27  for  three 
equal  Luttinger  parameters)  fixed  points  of  such  a  Y  junction 
in  the  (g\,g2,g3)  space,  and  compute  their  corresponding 
conductances  as  a  function  of  these  three  Luttinger  parameters. 
Moreover,  we  obtain  new  asymmetric  fixed  points,  which  are 
only  stable  for  nonidentical  TLLs,  thereby  providing  a  more 
complete  classification  of  the  conformally  invariant  BCs  for 
three  TLLs.  Such  asymmetric  fixed  points  have  been  identified 
using  perturbative  renormalization  group  analysis  in  the  weak 
coupling  regime.21,57,58 

Another  important  result  of  this  paper  is  a  direct  numerical 
verification,  through  density-matrix  renormalization  group 
(DMRG)  computations,65  of  the  analytical  predictions  for  the 
conductance  of  a  junction  of  two  nonidentical  wires.  Using  a 
recently  developed  method,66,67  which  allows  us  to  extract  the 
conductance  from  a  ground- state  calculation  in  a  finite  system, 
we  compute,  in  a  microscopic  lattice  model,  the  conductance  of 
a  junction  of  two  chains  of  interacting  spinless  electrons  with 
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FIG.  1.  (Color  online)  Junctions  of  two  and  three  TLLs  with 
different  Luttinger  parameters. 

different  interaction  strengths.  Our  numerical  results  show 
excellent  agreement  with  the  DEBC  predictions. 

The  outline  of  this  paper  is  as  follows.  In  Sec.  II,  we  set 
up  the  notation  and  present  the  model  in  the  bosonization 
framework.  In  Sec.  Ill,  we  present  the  DEBC  analysis  of 
a  junction  of  two  spinless  wires  and  show  that  the  scaling 
behavior  of  such  junctions  is  governed  by  a  single  effective 
Luttinger  parameter  ge  =  2 /(gf1  +  g2_1).  For  wires  with  spin- 
1/2  electrons,  we  obtains  two  effective  Luttinger  parameters 
gce,s  =  2/(1/ g[,s  +  \/gc{ s),  corresponding  to  charge  and  spin 
sectors.  Section  IV  contains  the  numerical  DMRG  calculations 
of  the  conductance  of  a  junction  of  two  spinless  chains  with 
different  interaction  strengths.  In  Sec.  V,  we  summarize  our 
results  on  the  Y  junction  with  the  detailed  DEBC  analysis 
presented  in  Appendix.  In  Sec.  VI,  we  present  the  analysis  of 
the  conductance  renormalization  when  the  wires  are  contacted 
to  Fermi-liquid  leads.  Finally,  we  conclude  in  Sec.  VII. 


II.  GENERAL  SETUP 


In  this  section,  we  present  the  model  in  the  bosonization 
framework  and  set  up  the  notation.  In  the  low-energy  limit,  the 
wires  are  described  by  TLLs  with  the  Euclidean  action1  • 
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where  g-  and  are,  respectively,  the  Luttinger  parameter 
and  the  plasmon  velocity  of  wire  /.  Different  wires  can 
have  different  electron-electron  interactions  and  consequently 
different  Luttinger  parameters  g;.  The  boson  fields,  <pi  and  ft, 
have  the  following  equal-time  commutation  relation, 


[<Pi(x)90j(x')]  =  inSijSgnix'  -  x),  (2.2) 

so  the  conjugate  momenta  of  ft  fields  are  given  by  = 
(dx(Pi)/2n.  Let  us  also  define  left-  and  right-moving  boson 
fields  as 


(Pi  =  4>f-  +  (t>?\  0i  =  (t>[  -  (t>*.  (2.3) 


where  rii  are  anticommuting  Klein  factors,  which  ensure 
the  correct  fermionic  statistics.  The  Klein  factors  play  no 
role  in  our  analysis  and  are  hence  neglected  throughout  the 
paper.  It  is  convenient  to  define  complex  variables  z  =  t  +  ix 
and  z  =  t  —  ix  such  that  the  left-  and  right-moving  bosons, 
corresponding  to  current  flowing  toward  and  away  from  the 
junction,  are,  respectively,  functions  of  z  and  2  only.  The  chiral 
current  operators  can  then  be  written  as 

jiR=z/^dA'  J>=z/kdA'  (2-5) 

where  ft  =  (ft  —  idx)/2  and  ft  =  ( ft  +  /ft)/ 2.  The  total 
current  is  proportional  to  the  difference  between  the  right  and 
the  left  currents:  ft  =  vi(JR  —  ftL). 

To  analyze  the  junctions  of  quantum  wires  with  unequal 
Luttinger  parameters,  it  is  convenient  to  introduce  rescaled 
bosonic  fields 


Si  =  Oi/Vgi’  <Pi  =  Vgi<Pi’  (2-6) 

which  effectively  have  a  noninteracting  action  [i.e.,  g  =  1  in 
Eq.  (2.1)].  Note  that  these  rescaled  fields  satisfy  the  original 
commutation  relations.  Similarly,  we  define  the  following 
rescaled  left-  and  right-moving  bosonic  fields: 

&•=#'  +  #';  0i=ff-4>?,  (2.7) 

in  terms  of  which  the  left-  and  right-moving  fermions  become 


ti 


L,R 


exp 


(2.8) 


In  the  absence  of  a  junction  (boundary),  the  correlation 
functions  of  the  rescaled  fields  are  given  by 


8ij 

(<Pi(z,z)<Pj(w,w))  =  — |r  ln[(z  -  w)(z  -  w)], 

(- §i(z,z)0j(w,w ))  =  -  Y  In [(z  -  w)(z  -  U))L  (2.9) 
{(pi(z,z)0j(w,w))  =  0. 


Imaginary  time  ordering  is  implied  for  all  correlation  functions 
here  and  throughout  the  paper.  Different  bosonic  fields  above 
are  of  course  uncorrelated.  However,  as  we  will  show  below, 
the  presence  of  a  junction  mixes  these  fields  and  effectively 
reduces  the  independent  bosonic  degrees  of  freedom  by  half. 

At  the  end  of  wires,  x  =  0,  fermions  can  hop  between 
different  wires,  this  process  is  described  by  a  single  particle 
hopping  Hamiltonian 


Hb  =  -  J2[tijeia‘^j(0)fj(0)  +  H.c.],  (2.10) 

ij 


The  fermions  ift(x)  =  elkFX\j/lR(x)  +  e~lkpx\l/R(x),  with 
iftL(x)  and  ^sR(x)  the  linearized  left-  and  right-moving 
fermionic  fields,  can  be  written  in  terms  of  the  above  bosons 
through 


JfL,R  =  m  _  Hi  ci(<Pi±0i)/V2’ 

1  \f2n  \f2n 


(2.4) 


where  ft  and  eft  are  the  strength  and  the  phase  of  the 
hopping  amplitude  between  wires  /  and  j.  Without  loss  of 
generality22,27  (at  least  for  junctions  of  two  or  three  wires), 
we  only  consider  the  symmetric  case  ft  =  /  in  this  paper. 
The  phases  encode  the  distribution  of  magnetic  fluxes  at  the 
junction,  which  for  a  junction  of  three  or  more  wires  can  play 
a  crucial  role  in  the  properties  of  some  RG  fixed  points.22,27 
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III.  JUNCTION  OF  TWO  WIRES:  DELAYED  EVALUATION 
OF  BOUNDARY  CONDITIONS  ANALYSIS 

In  this  section,  we  analyze  the  stability  of  RG  fixed  points, 
and  compute  their  corresponding  conductances  for  a  junction 
of  two  TLL  quantum  wires  with  unequal  Luttinger  parameters. 
By  using  the  DEBC  method  of  Ref.  27,  we  show  that  the 
properties  of  such  junctions  only  depend  on  an  effective 
Luttinger  parameter 

Sp  =  {gP+gP)/ 2-  (3.1) 

Briefly,  the  junction  of  two  wires  is  totally  decoupled  for 
ge  <  1,  and  has  a  conductance  ge(e2/h)  for  ge  >  1.  This 
result  is  consistent  with  what  was  found  in  Ref.  43,  where  the 
tunneling  between  fractional  quantum  Hall  edge  states  with 
different  filling  fractions  was  discussed.  The  DEBC  method 
used  in  this  paper,  however,  has  the  advantage  that  it  can  be 
straightforwardly  generalized  to  junctions  of  more  than  two 
wires  (see  Sec.  V  and  Appendix). 


A.  Delayed  evaluation  of  boundary  conditions  method 

For  one-dimensional  quantum  impurity  problems,  one  often 
invokes  the  conformal  symmetry  of  the  bulk  system  (Luttinger 
liquids  in  our  case)  and  assumes  that  the  effect  of  an  impurity 
(junction),  at  low  energies,  is  imposing  a  conformally  invariant 
boundary  condition  (BC),  which  describes  the  renormalization 
group  (RG)  fixed  point.  This  methodology  of  relating  the  BC 
and  RG  fixed  points  of  the  system  is  called  boundary  conformal 
field  theory  (BCFT),  and  has  proved  greatly  successful  in  the 
study  of  quantum  impurity  problems.63,64 

A  useful  technique  within  the  framework  of  BCFT  is  the 
DEBC  method,  which  hugely  simplifies  the  evaluation  of  the 
scaling  dimensions,  Aqb,  of  boundary  operators,  Ob,  with  a 
given  BC.27  The  scaling  dimension,  in  turn,  determines  the 
leading  scaling  behavior  of  a  given  operator  under  the  RG 
flow,  and  thus  govern  the  stability  of  the  RG  fixed  points: 
In  general,  an  RG  fixed  point  (boundary  condition)  is  stable 
if  all  boundary  operators  are  either  equivalent  to  identity  or 
irrelevant  A oB  >  1.  Moreover,  the  conductance  associated 
with  the  given  fixed  point  can  be  readily  computed  from 
the  BC.  For  a  complete  description  of  the  DEBC  method, 
we  refer  the  reader  to  Refs.  27  and  29.  Here,  we  simply 
apply  this  method  to  a  junction  of  two  quantum  wires  with 
unequal  Luttinger  parameters,  commenting  only  on  some  key 
ingredients. 


B.  Junctions  of  two  quantum  wires 

The  first  step  of  the  DEBC  method  is  to  write  an  ansatz 
for  the  conformally  invariant  BCs  describing  the  RG  fixed 
points.  The  next  step  is  to  list  all  the  boundary  operators, 
which  can  possibly  become  relevant  and  make  the  fixed  point 
unstable,  and  compute  their  scaling  dimensions  with  such 
ansatz  for  every  point  in  the  parameter  space,  g\  and  g2-  If 
for  a  given  ansatz,  none  of  these  boundary  operators  have  a 
scaling  dimension  smaller  than  one  (in  some  region  of  the 
parameter  space  known  as  the  stability  region),  we  have  found 
a  stable  RG  fixed  point.  In  case  of  a  junction,  a  natural  ansatz 
can  be  expressed  in  terms  of  a  rotation  matrix  7 Z  that  relates 


outgoing  to  incoming  bosonic  fields 

$R  =  HQ1,  (3.2) 

where  <j>LR  =  (4>i'r,  . . .  ,<pR’R)T  are  /-component  vector 
fields. 

The  most  important  boundary  operators,  in  case  of  a 
junction  of  two  wires,  correspond  to  the  following  processes: 
tunneling  of  chiral  fermions  between  the  two  wires,  and 
backscattering  within  the  individual  wires.  It  is  useful  to  in¬ 
troduce  a  compact  notation  for  boundary  operators  describing 


the  single-particle  tunneling  processes, 

6 

II 

H 

III 

a 

(3.3) 

where  a,b  =  R,L. 

We  can  then  list  these  six  fundamental 

boundary  operators 

in  terms  of  the  rescaled  boson  fields  as 

follows: 

T*  RL 

1 21(12) 

.  e±^^)eT2^+^\ 

’-pRL  0 
1 11(22) 

(3.4) 

tLL(RR) 

1 21 

The  boundary  operators  corresponding  to  multiparticle  pro¬ 
cesses  are  not  forbidden  and  can  be  generated  as  higher-order 
perturbation  processes  even  they  are  not  presence  in  the 
bare  Hamiltonian.  In  general,  they  can  be  constructed  from 
these  fundamental  boundary  operators  and  have  larger  scaling 
dimensions  and  are  less  relevant  than  the  single-particle 
processes. 

All  the  above  boundary  operators  have  the  generic  form 
Ob  ~  eia-v+lb'0 ?  where  a,b  are  vectors  that  contain  the 
prefactors  of  the  (pi  and  0t  fields.  By  eliminating  the  redundant 
degrees  of  freedom  with  Eq.  (3.2),  and  using  Eq.  (2.9),  the 
scaling  dimension  of  the  generic  Ob  above  can  then  be  written 
in  terms  of  7 Z  as 


A%B  =  \\KT(a-b)  +  (a  +  b)\2,  (3.5) 

where  the  superscript  T  represents  matrix  transpose. 

To  find  all  the  7 Z  matrices  that  correspond  to  stable  fixed 
points,  it  is  convenient  to  express  (pi  fields  in  terms  of  the 
following  zb  fields: 


<P+ 


+  \/g202  m 


(p- 


Vg2<Pl  ~  Vglfc 


(3.6) 


y/gl  +  §2  y/gl  +  g2 

Corresponding  6±  are  defined  in  a  similar  manner.  The  six 
fundamental  boundary  operators  can  then  be  written  as 


’-pRL 

12\{\2) 

rLL{RR ) 
1 21 


±i^= 

e  , 


’-pRL 
1 11(22) 


f. J  *fge  e±iyf8~e0-  . 


e±i^fg~e0-  , 


(3.7) 


Here,  we  have  dropped  all  el6+  terms  as  they  are  effectively 
an  identity  at  the  boundary:  Charge  conservation  requires 
J]-  JfR  —  JtL  =  0,  which,  using  Eqs.  (2.5)  and  (2.6),  gives 
3 T  0+  =  0  (i.e.,  Dirichlet  BC  on  0+),  and  makes  el0+  an  effective 
identity.27  The  simplified  boundary  operators  in  Eq.  (3.7)  then 
only  depend  the  effective  Luttinger-liquid  parameter  ge. 

In  terms  of  the  left-  and  right-moving  =b  fields  defined  in  a 
similar  manner  to  Eq.  (2.7),  the  Dirichlet  BC  on  0 +  gives 


0 


R 

+ 


L  I 

+  lx=0" 


(3.8) 


075451-3 


HOU,  RAHMANI,  FEIGUIN,  AND  CHAMON 


PHYSICAL  REVIEW  B  86,  075451  (2012) 


Now  we  only  need  to  specify  the  BC  relating  the  fields 
as 


0^  =  (3.9) 

Because  there  is  a  single  pair  of  RL  fields,  only  the  Neumann- 
BC  (N-BC)  and  the  Dirichlet-BC  (D-BC),  K*L'D  =  ±1,  are 
allowed.  By  using  Eq.  (3.5)  on  boundary  operators  listed  in 
Eq.  (3.7),  we  obtain  the  scaling  dimension  of  each  operator 
with  the  N-BC  and  D-BC. 


Ob 

a£b  (N-BC) 

Ags  (D-BC) 

n rRL  ’-rRL 

12\  ’  1 12 

1  /ge 

0 

n rRL  ’-rRL 

1  11  »  122 

0 

8e 

n rLL  ’-rRR 

221  >  12\ 

1  /ge 

ge 

Here,  the  scaling  dimension  0  indicates  that  Ob  is  equivalent 
to  identity  operator  1  for  the  given  boundary  condition. 

From  the  table  above,  we  conclude  that  the  N-BC  is  stable 
when  ge  <  1  and  D-BC  is  stable  when  ge  >  1.  As  shown 
below,  the  N-BC  corresponds  to  a  fixed  point  where  two  wires 
are  disconnected  and  the  D-BC  corresponds  to  a  fixed  point 
where  two  wires  are  maximally  connected  with  conductance 
Gd  =  ge(e2 /h).  Let  us  recall  the  well-known  results  of  the 
Kane  and  Fisher  problem,  namely,  a  single  impurity  in  a 
spinless  TLL  wire,  which  is  equivalent  to  a  junction  of  two 
wires  with  the  same  Luttinger  parameter,  g.u  There,  two 
RG  fixed  points  are  identified,  a  totally  disconnected  fixed 
point  for  g  <  1  and  a  maximally  connected  fixed  point  for 
g  >  1  with  GD  =  g{e2  / h).  As  ge  -*  g  when  gi  =  g2  =  g, 
the  fixed  points,  which  we  identified,  generalize  Kane  and 
Fisher’s  results  to  the  case  of  two  quantum  wires  with  different 
Luttinger  parameters. 

A  similar  DEBC  analysis  can  also  be  applied  to  a  junction 
of  two  quantum  wires  with  spin- 1/2  electrons.29  In  that 
case,  we  obtain  two  effective  Luttinger  parameters  gce's  — 
2/(1  /g\S  +  l/#2  ' s)>  corresponding  to  charge  and  spin  sectors, 
which  govern  the  scaling  behavior  of  boundary  operators  near 
the  fixed  points.68  As  a  consequence,  the  stable  fixed  points 
of  such  system  are  connected  to  those  of  equal  Luttinger 
parameters  but  with  gc" s  — >  gce,s .1546,18 


C.  Conductance  for  each  fixed  point 

Here,  we  will  compute  the  conductance  associated  with 
N-BC  and  D-BC.  The  conductance  tensor  is  defined  through 
the  following  current- voltage  relation: 

4=EG‘VV’  (3-10) 

j 

where  current  is  defined  as  positive  when  flowing  toward 
junction.  In  the  linear-response  regime,  the  conductance  above 
can  be  evaluated  via  the  Kubo  formula,27 

e2  1  f°°  rL 

Gij  =  lim  - - /  dzelcoT  /  dx(Ji(y,T)Jj(x,0))- 

o)^0+  h  coL  J_00  Jo 

(3.11) 


The  current  correlation  function  can  be  rewritten  as  a  sum  of 
chiral-current  correlation  functions  as 

0)}  =  (jR{y,T)jf{xfi))  +  (j^y,T)jf{x,Q)) 

-  r)jj(x,0))  -  {jj(y,x)jf(x,0% 

(3.12) 


where  the  left  and  right  currents  are  defined  in  Eq.  (2.5).  By 
using  Eq.  (2.6),  we  can  write  the  chiral  currents  in  terms  of 
rescaled  boson  fields  as  follows: 


JiR  =  +i^d&  =  -  V*  j;*, 

J f  =  =  vW". 


(3.13) 


where  we  have  defined  left  and  right  currents  associated  with 
the  rescaled  fields. 

To  evaluate  the  above  correlation  functions  given  a  BC,  it  is 
convenient  to  first  express  the  boundary  conditions  1Z  directly 
in  the  rescale  boson  field,  4>^,R,  basis.  From  Eq.  (3.8)  and 
Eq.  (3.9),  we  can  derive  the  rotation  matrix 


nN  =  i;  nD 


(g\-gi 
g\+g2 

g\+gi 


g\+gi 
g2  g  1 

g\+g2 


(3.14) 


for  the  N-BC  and  D-BC,  respectively.  From  Eq.  (3.13),  the 
left  and  right  currents  are  constrained  by  the  JR( 0,r)  = 
TZij  BC  at  the  origin,  which,  upon  unfolding  the 

current,  implies 

JiR(x,T)  =  nijg(-x,r),  (3.15) 

for  x  >0.  Therefore,  all  the  right-moving  currents  can  be 
interpreted  as  left-moving  on  the  x  <  0  domain. 

Now,  the  chiral  current  correlation  functions  can  be  evalu¬ 
ated, 


(JfiMJfizj)) 


(jR{Zi)jf(Zj)) 


(jf{Zi)jf{Zj)) 


Sij 

gi 

An2 

(■ Zi 

-Zj)2 

Sij 

gi 

An2 

C zt 

-Zj)2 

IZij 

\/gigj 

An2 

(■ zi 

-Zj)2 

TZji 

VsisJ 

An2 

(■ Zi 

-Zj)2 

(3.16) 


By  inserting  these  correlation  functions  into  the  Kubo  formula 
Eq.  (3.11),  and  after  some  algebra,  we  obtain  a  concise  relation 
between  conductances  and  boundary  conditions27 


Gij  —  ^  V SiSj  ($ij  ij )• 


(3.17) 


With  the  N-BC  and  D-BC  represented  by  the  rotation  matrices 
of  Eq.  (3.14),  we  immediately  conclude  that 


(3.18) 


As  expected,  the  N-BC  corresponds  to  a  fixed  point  with 
decoupled  wires  and  D-BC  corresponds  to  a  fixed  point  with 
conductance  ge(e2/h). 
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Here,  we  shall  emphasize  that  the  correlation  functions 
listed  in  Eq.  (3.16)  include  only  the  universal  part  for  a  given 
boundary  condition.  There  are  also  nonuniversal  contributions 
to  the  correlation  functions,  which,  in  general,  decay  faster  and 
become  irrelevant  at  long  distance.  However,  when  the  univer¬ 
sal  part  vanishes,  which  is  the  case  for  correlations  between 
different  wires  with  the  N-BC,  the  higher-order  contributions 
dominate  and  could  lead  to  nonlinear  conductance.  In  the 
next  section,  we  will  perform  the  DMRG  calculations  and 
confirm  that  Eq.  (3.16)  indeed  represents  the  universal  part  of 
chiral-current  correlation  functions. 


IV.  JUNCTION  OF  TWO  WIRES:  DENSITY-MATRIX 
RENORMALIZATION  GROUP  CALCULATIONS 

In  this  section,  we  perform  numerical  computations  of 
the  conductance  of  a  junction  of  two  Luttinger  liquids  with 
different  Luttinger  parameters  g\  and  A  microscopic  lattice 
model  with  the  Luttinger-liquid  physics,  which  is  suitable  for 
numerical  calculations,  is  the  one-dimensional  tight-binding 
model  of  interacting  spinless  electrons, 

Hi  =  ^  '  ^i,mCi,m+ 1  4“  H.C.  +  V(  (fli  m  ~ 

m  ' 

(4.1) 

for  wire  /,  where  the  hopping  amplitude  is  set  to  unity.  A 
junction  of  two  wires  can  be  described  by 


We  extract  the  conductance  Gtj  from  the  following  asymp¬ 
totic  (large  x)  relationship  in  the  finite  system  shown  in 
Fig.  2(b): 


(■ J‘r(x)Jjl(x )) 


41  sin 


(4.3) 


where  (JlR(x)j[(x))  is  the  ground- state  correlation  function  of 
chiral  currents  in  wires  i  and  j  in  Fig.  2(b).  In  terms  of  total 
charge  N  and  current  J  operators,  which  can  be  modeled  on 
the  lattice,  we  generically  have67 

(4(x)//(x))  =  (x)jj(x))  -  ^~{N\x)Jj{x)). 

'  '  2  vtVj  2  vj 

(4.4) 


Note  that  in  the  time-reversal  symmetric  case  considered  here, 
the  second  term  in  the  expression  above  vanishes,  and  we  only 
need  to  compute  a  static  current-current  correlation  function 
(Jl(x)Jj  (x)).  In  terms  of  the  lattice  creation  and  annihilation 
operators  appearing  in  Hamiltonian  (4.1),  we  have 

J1  (m  +  —  c\mCi^m+ 1). 

All  we  need  to  do  now  is  to  numerically  compute  {J1  (x)J2(x)) 
for  the  above  current  operator  and  divide  it  by  2v\ u2  to 
obtain  (JR(x)JR(x)).  The  numerical  calculations  are  done  for 
a  system  of  180  sites  in  each  of  the  two  wires.  The  truncated 
number  of  states  in  our  DMRG  computations  is  1100. 


H  =  Hi  +  H2  -  tc\  0c2, o  -  f<4,oci.o-  (4-2) 

The  parameters  gt  and  of  the  Luttinger-liquid  Hamiltonian 
(2.1)  are  related  to  the  interaction  strength  V  through  the  Bethe 
ansatz  (see,  e.g.,  Refs.  66  and  69).  We  use  the  method  of 
Refs.  66  and  67  to  compute  the  conductance.  This  method 
allows  us  to  extract  the  conductance  from  a  ground- state  static 
calculation  in  a  finite  system  as  explained  below.  The  semi¬ 
infinite  junction  of  Hamiltonian  (4.2)  is  depicted  in  Fig.  2(a), 
while  the  corresponding  finite  system  used  in  the  numerics  is 
shown  in  Fig.  2(b). 
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FIG.  2.  (Color  online)  (a)  The  junction  of  two  semi-infinite 
Luttinger-liquid  wires  with  interaction  strength  V]  and  V2.  Hopping 
amplitude  is  set  to  unity  in  the  bulk  of  the  the  two  wires,  and  is  equal 
to  t  at  the  junction,  (b)  The  corresponding  finite  system  used  for 
DMRG  calculations  according  to  the  method  of  Refs.  66  and  67. 


A.  Repulsive  effective  interaction  ge  <  1 

For  ge  <  1,  we  have  Gu  =  0,  which  implies  that  the 
leading  term  [Eq.  (4.3)]  in  the  (JR(x)JR(x))  correlation 
function  vanishes.  If,  as  a  function  of  ^  sin(|x),  the  computed 
(JR(x)Jl(x))  decays  faster  than  a  power  law  with  exponent 
—2,  we  have  a  signature  of  a  vanishing  G12.  Our  numerical 
results  indeed  confirm  this:  for  any  combination  of  g\ 
and  g2  with  ge  <  1,  we  find  that  (JR(x)Jl(x))  decays  as 
^  sin(jx)~a(^8e>)  with  a(ge)  >  2.  The  exponent  a  only  depends 
on  ge ,  and  is  independent  of  individual  gt  and  the  hopping 
amplitude  t.  The  prefactor  of  the  correlation  function  depends 
on  the  hopping  amplitude  t  since  universality  is  a  property 
of  the  leading  term,  and  the  coefficient  of  the  subleading 
term  observed  here  can  depend  on  microscopic  details  such 
as  t. 

Note  that  the  correlation  functions  for  different  combina¬ 
tions  of  g\  and  g2,  which  have  the  same  ge ,  collapse  not  only 
in  the  large  x  limit  but  also  close  to  the  microscopic  length 
scales.  This  strongly  suggests  that  the  single  parameter  ge 
determines  all  the  subleading  corrections  to  the  correlation 
function.  This  behavior  can  be  understood  by  noting  that 
all  the  boundary  operators  in  Eq.  (3.7)  depend  on  ge  (as 
opposed  to  individual  g\  and  g2)  after  dropping  the  (p+  and 
6+  fields  due  to  current  conservation.  Hence,  the  measured 
correlation  functions  should  also  be  determined  only  by  the 
effective  Luttinger  parameter  ge  and  the  hopping  strength  t. 
Interestingly,  the  exponent  a  is  indeed  close  to  the  scaling 
dimension  of  the  leading  irrelevant  operator  (i.e.,  2 / ge)  but  we 
are  not  able  to  make  a  definitive  statement  due  to  finite-size 
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£  =  180,  ge  =  0.83 


(b)  £  =  180,  ge  =  0.87 


FIG.  3.  (Color  online)  The  correlation  function  (J^(x)Jl(x))  for 
different  combinations  of  g i  and  g2  corresponding  to  two  values 
of  ge  <  1.  The  data  exhibits  very  good  collapse,  and  indicates  a 
vanishing  conductance. 


effects  and  limited  numerical  precision.  It  is  worth  mentioning 
that  we  have  also  considered  combinations  of  g\  and  g2  where 
one  wire  has  attractive  interactions  (gi  >  1),  but  this  does  not 
affect  the  behavior  of  the  junction  as  long  as  ge  <  1. 

The  results  are  summarized  in  Fig.  3.  We  have  considered 
two  values  of  ge  =  0.871,0.83  and  four  combinations  of  g\ 
and  g2  for  each  ge  as  shown  below. 


Vi 

v2 

gi 

g2 

ge 

0.463 

0.463 

0.871 

0.871 

0.871 

0 

0.9 

1 

0.771 

0.871 

0.347 

0.576 

0.9 

0.843 

0.871 

-0.285 

1.145 

1.1 

0.720 

0.871 

0.632 

0.632 

0.830 

0.830 

0.830 

0 

1.2 

1 

0.709 

0.830 

0.347 

0.904 

0.9 

0.770 

0.830 

-0.285 

1.415 

1.1 

0.666 

0.830 

Here,  we  have  considered  two  different  values  of  hopping 
amplitude,  t  =  0.5, 0.7. 

B.  Attractive  effective  interaction  ge  >  1 

In  this  case,  we  expect  any  tunneling  amplitude  t  ^  1  at 
the  junction  to  heal  and  result  in  a  universal  conductance 
gee2  /  h.  In  our  numerics,  we  have  used  two  values  of  tunneling 
amplitude,  t  =  0.7, 0.9.  Similarly  to  the  repulsive  case,  we 
consider  two  values  of  ge  =  1.175,1.258;  for  each  ge,  we 
consider  four  combinations  of  g\  and  g2  shown  below. 
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g2 

ge 

-0.4625 

-0.4625 

1.175 

1.175 

1.175 

0 

-0.9 

1 

1.423 

1.175 

0.347 

-1.19 

0.9 

1.690 

1.175 

-0.285 

-0.637 

1.1 

1.260 

1.175 

-0.632 

-0.632 

1.258 

1.258 

1.258 

0 

-1.2 

1 

1.690 

1.258 

0.347 

-1.460 

0.9 

2.087 

1.258 

-0.285 

-0.960 

1.1 

1.468 

1.258 

Again,  we  obtain  very  good  agreement  with  the  analytical 
predictions.  The  results  are  summarized  in  Fig.  4.  With  fixed 
ge,  the  data  points  collapse  for  different  combinations  of  g\ 
and  g2  and  different  values  of  hopping  t.  The  numerical  results 
approach  the  analytical  prediction  for  a  conductance  of  gee2  /  h 
(solid  black  line)  in  the  asymptotic  limit. 

V.  JUNCTION  OF  THREE  WIRES 

A  junction  of  three  quantum  wires  with  equal  Luttinger 
parameters  has  three  distinct  types  of  fixed  points  described 
by  a  rotation  matrix  ansatz:  decoupled  fixed  point,  chiral- x± 
fixed  points,  and  Dirichlet  fixed  points,  which  are  respectively 
stable  forg<l,l<g<3  and  g  >  3.  (There  is  an  additional 
less  understood  time-reversal-invariant  M  fixed  point  for 
1  <  g  <  3,  which  we  do  not  consider  in  this  work.)22,27,66,67 
Thus,  it  is  reasonable  to  expect  that  these  fixed  points  remain 
stable  if  we  slightly  change  the  Luttinger  parameters  and  make 
them  unequal.  Here,  we  first  determine  this  region  of  stability 
around  the  g\  =  g2  =  g3  line.  In  addition,  we  identify  three 


£=  180,  ge  =  1.175 


(b)  1=  180,  ge  =  1.258 


ln[fsm(?x)] 

FIG.  4.  (Color  online)  The  correlation  function  (Jr(x)JI(x))  for 
different  combinations  of  g i  and  g2  corresponding  to  two  values 
of  ge  >  1.  The  data  exhibits  very  good  collapse,  and  indicates  a 
conductance  of  gee2/h  because  of  the  asymptotic  agreement  with  the 
exact  theoretical  prediction  shown  by  a  sold  black  line. 
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TABLE  I.  The  scaling  dimensions  of  leading  irrelevant  boundary  operators  and  the  conductance  tensor  for  each  stable  fixed  point  of  a  Y 
junction.  The  detailed  analysis  for  obtaining  these  scaling  dimensions  is  given  in  Appendix  with  the  corresponding  operators  listed  in  Table  II. 
The  asymmetric  fixed  point  At  represents  a  boundary  condition  where  the  wire  i,  for  i  =  1,2,3,  is  decoupled  from  the  junction.  Here,  the 
scaling  dimensions  of  all  leading  irrelevant  operators  run  over  the  indices  for  all  possible  combinations  of  i  ^  j  ^  k.  We  have  also  introduced 
following  notations:  the  cyclic  identification  g0  =  g3  and  g4  =  gu  the  two  indices  antisymmetric  tensor  £jj±\  =  ±1  and  0  otherwise;  and  the 
index  m  satisfies  m  ^  j  ^  k.  The  conductance  tensors  are  given  units  of  e2/h  and  are  defined  through  Ij  =  GjkVk. 


Fixed  point 

Scaling  dimensions  A 

Conductance  Gjk  [e2/h\ 

Decoupled  (N-BC) 

(gi  +  g.i)/2gigj 

0 

Dirichlet  (D-BC) 

giigj  +  gk)/2(gi  +  gl  +  g3) 

2  [gjSjk  —  gjgk/(gi  +  g2  +  £3)] 

Chiral  (x±-BC) 

2  giigj  +  gk)/(glglg3  +  gl  +  gl  +  g3 ) 

O  8j (SI  +<?2+£3)<S jk+g j 8k(T8m ejk- 1) 

*1  +87+fn 

(gl  +  g2  +  g3  +  gig2g3)/2gi(gi+\  +  gi- 1) 

^^(-1  +  %  +  8,k  +  2 8Jt  -  3 

Asymmetric  (ArBC) 

2gi+igi-\/(gi+i  +gi- 1) 

2(gl  +  g2  +  g3)/gi(gi+l  +  gi- 1) 

asymmetric  fixed  points,  only  realized  for  unequal  Luttinger 
parameters,  in  which  one  of  the  wires  is  decoupled  from  the 
junction  and  the  other  two  wires  are  fully  connected.  These 
asymmetric  fixed  points  have  important  consequences  for  the 
stability  of  N,  D,  and  x  ones:  there  are  regions  of  the  parameter 
space  near  the  transition  points  gi  =  1,3  on  the  g\  =  g2  =  g3 
line  where  small  perturbations  normal  to  the  equal-g  line 
would  drive  the  system  into  one  of  the  asymmetric  fixed 
points. 

In  Table  I,  we  first  summarize  the  scaling  dimensions  of  the 
leading  boundary  operators  as  well  as  their  corresponding  con¬ 
ductances  as  a  function  of  the  three  Luttinger  parameters  for 
each  boundary  condition.  These  scaling  dimensions  determine 
the  stability  of  the  fixed  point:  when  all  them  are  larger  than 
one  in  certain  parameter  region,  the  given  fixed  point  is  stable. 
From  Table  I,  we  observe  that  the  three  scaling  dimensions  of 
the  leading  irrelevant  operators  for  the  asymmetric  fixed  points 
are  each  the  inverse  of  those  of  the  decoupled,  Dirichlet  and 
chiral  fixed  points.  Hence,  in  any  given  point  of  the  (gi  ,gi,g?>) 
parameter  space,  there  exists  at  least  one  stable  fixed  point.  In 
other  words,  the  decoupled,  Dirichlet,  chiral,  and  asymmetric 
fixed  points  fully  cover  the  phase  diagram  of  a  Y  junction  of 
spinless  TLL  wires. 

In  this  paper,  the  DMRG  analysis  is  not  applied  to  the 
junction  of  three  quantum  wires  with  unequal  Luttinger 
parameters.  Such  analysis  with  equal  Luttinger  parameters 
was  performed  in  Ref.  67,  but  is  beyond  the  scope  of  the 
present  paper.  We  discuss  the  physical  properties  of  each 
of  these  stable  fixed  points  in  the  remainder  of  this  section. 
The  detailed  analysis,  based  on  the  DEBC  method,  can  be 
found  in  Appendix.  As  the  conductances  of  fixed  points  can 
be  evaluated  in  the  similar  way  as  in  Sec.  IIIC,  we  simply 
write  down  the  results  in  the  following  discussion. 

A.  Decoupled  fixed  point 

The  decoupled  fixed  point  corresponds  to  the  Neumann 
BC  for  all  the  bosonic  fields  (p .  From  Table  I,  we  see  that 
the  scaling  dimensions  of  the  leading  irrelevant  operators  are 
equal  to  l/gle,J  for  gle,J  =  2gigj/(gi  +  gj),  which  is  the  same 
as  Eq.  (3.1)  for  a  pair  of  wires.  Hence,  the  decoupled  fixed 
point  (N-BC)  is  then  stable  when  the  N  fixed  point  is  stable 
for  all  three  possible  pairs  of  wires.  One  can  simply  check  that 


these  scaling  dimensions  reduce  to  l/g  when  the  Luttinger 
parameters  are  all  equal.  In  Fig.  5,  the  stability  region  of  the 
decoupled  fixed  point,  A@b  >  1,  is  shown  in  red.  As  expected, 
the  conductance  of  the  decoupled  fixed  points  is  simply 

=  0.  (5.1) 


B.  Dirichlet  fixed  point 

As  discussed  in  Appendix,  one  can  construct  three  inde¬ 
pendent  linear  combinations  of  the  bosonic  field  (p  such  that, 
akin  to  (p+  field  for  junctions  of  two  wires,  one  of  them,  known 
as  the  center  of  mass  field,  always  satisfies  the  Neumann  BC 
due  to  charge  conservation.  The  Dirichlet  fixed  point  corre¬ 
sponds  to  imposing  the  D-BC  on  the  other  two  combinations. 
None  of  single-particle  processes  becomes  identity  with  such 
boundary  condition.  Instead,  some  of  two-  or  more-particle 
processes  reduce  to  identity,  which  suggests  that  the  Dirichlet 
fixed  point  is  associated  with  a  certain  type  of  Andreev 
reflection  that  enhances  the  conductance. 

From  Table  I,  the  scaling  dimensions  of  all  leading 
irrelevant  operators,  Agc  =  g,(gj  +  gk)/^(g\  +  gi  +  gs)  V 
i  z/z  j  z/z  k ,  reduce  to  g/ 3  when  g\  =  g2  =  #3  =  g-  Hence, 
the  D-BC  becomes  stable  at  g  >  3,  consistent  with  Ref.  27. 
By  requiring  Aqb  >  1,  we  obtain  the  stability  region,  painted 
in  green  in  Fig.  5,  of  the  Dirichlet  fixed  point  for  unequal 
Luttinger  parameters. 

The  conductance  of  the  Dirichlet  fixed  point  is  given  by 


G 


D 

jk 


8jSjk  ~ 


8j_gk 

gl  +  g2  +  #3  _ 


(5.2) 


where  we  have  made  a  cyclic  identification  g0  =  #3  and  g4  = 
gi.  When  all  the  Luttinger  parameters  are  equal,  we  have 
Gfk  =  g(e2 /  h)(28jk  —  2/3),  which  reproduces  the  result  in 
Ref.  27. 


C.  Chiral-x±  fixed  points 

The  chiral- x±  fixed  points  have  a  particular  transport 
feature:  the  realization  of  x+  or  X-  fixed  points,  with  the 
incoming  current  respectively  flowing  clockwise  or  counter¬ 
clockwise  into  one  of  the  adjacent  wires,  depends  on  the 
direction  of  the  threaded  magnetic  field  into  the  ring.27  (This 
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g3t  8  =  (gi  +g2  +  g3>/3 


g3  =  0 


3V2 g 


The  conductances  for  chiral- x±  fixed  points  are,  in  turn, 
give  by 


G 


x± 

jk 


8j(8 1  +  82  +  83)^ jk  +  8j8k(rF8m€ jk  —  1) 
h  glg2g3  +  gl  +  82  +  83 


(5.3) 


(b)  .  g  =  V3/2 

(c)  A  g  =  1 

mL 

(e)  A  1  =  3 

fh. 

(f)  A  8  = 4 

A/4A\ 

FIG.  5.  (Color  online)  We  plot  the  stability  regions  of  the  different 
fixed  points  in  the  triangular  cross  sections,  shown  in  upper  panel,  of 
fixed  g  =  ( g\  +  g2  +  g3)/3  in  Luttinger  parameters  space,  g i,  g2 ,  and 
g3.  In  all  panels,  the  decoupled,  chiral,  and  Dirichlet  fixed  points  are 
painted,  respectively,  in  red,  orange,  and  green,  and,  the  asymmetric 
fixed  points  Ai,2j3  are  pained  in  yellow,  gray,  and  blue,  respectively. 
The  white  areas  represent  an  overlap  of  two  asymmetric  fixed  points 
and  the  equal-g  points  are  indicated  by  a  black  dot  at  the  center  of 
triangle,  (a)  For  g  <  2/3,  only  the  decoupled  fixed  point  is  realized. 

(b)  The  x±  and  Dirichlet  fixed  points  are  not  stable  when  g  <  1. 

(c)  The  point  g i,2,3  =  1  is  an  exactly  marginal  point  surrounded 
by  asymmetric  fixed  points,  (d)  The  x±  fixed  points  appear  at  the 
center  of  cross  section  for  1  <  g  <  3,  surrounded  by  asymmetric 
fixed  points.  The  stability  region  for  the  N-BC  are  pushed  to  the 
corners  where  two  of  the  three  Luttinger  parameters  become  much 
less  than  one.  (e)  The  point  £1,2,3  =  3  is  another  exactly  marginal 
point,  again  surrounded  by  asymmetric  fixed  points.  The  chiral  fixed 
points  have  extremely  small  stability  regions  (difficult  to  see  in  the 
figure)  located  between  any  two  asymmetric  fixed  points,  (f)  The 
Dirichlet  fixed  point  appears  at  the  center  of  triangle  for  g  >  3.  Note 
that  the  figure  is  not  a  schematic,  and  represents  the  exact  domain 
boundaries  for  the  shown  g. 


point  will  become  more  apparent  in  the  next  section  when 
discussing  Fermi-liquid  leads)  When  g\  =  g2  =  g3  =  g,  the 
scaling  dimensions  listed  in  Table  I  for  both  x±  fixed  points 
reduce  to  ^ ,  and  hence  x±  fixed  point  is  stable  for  1  <  g  <  3, 
which  is  consistent  with  what  found  in  Ref.  27.  Again,  we 
obtain  the  stability  region  of  the  chiral  fixed  points  through 
>  1.  In  Fig.  5,  such  region  is  shown  in  orange. 


where  €jJ±l  =  ±1  while  ejk  =  0  for  j  =  k ,  and  gm  /  gj,gk. 


D.  Asymmetric  fixed  points 

The  asymmetric- Az  fixed  points  have  the  property  that  wire 
i  is  decoupled  from  the  junction.  Such  fixed  points  are  a  new 
feature  of  a  system  with  unequal  Luttinger  parameters.  As 
shown  in  Table  I  and  mentioned  earlier,  the  scaling  dimensions 
of  the  leading  irrelevant  operators  of  A/  fixed  points  are 
inverse  to  those  of  the  decoupled,  Dirichlet  and  chiral  fixed 
points.  Hence,  the  stability  regions  of  A/  fixed  points  are 
complimentary  to  those  of  the  other  fixed  points,  and  the  entire 
parameter  space  is  covered  by  at  least  one  of  the  stable  fixed 
point  presented  in  this  work.  In  Fig.  5,  the  stability  regions 
of  A  1,2,3  fixed  point  are  shown  in  yellow,  gray,  and  blue, 
respectively.  The  regions  where  two  asymmetric  fixed  points 
overlap  are  shown  in  white. 

The  conductances  of  the  asymmetric  fixed  points  A/  are 
give  by 

G%  =  g‘e+u~l  j(~l  +  Su  +  8ik  +  2 Sjk  -  3 SijSik),  (5.4) 

where  g™,n  is  the  effective  Luttinger  parameter  for  the  pair  of 
wires  m  and  n.  To  give  an  idea  of  the  properties  of  asymmetric 
fixed  points,  we  hereby  write  the  conductance  tensor  explicitly 
for  Ai, 


/O 

0 

0  ) 

CA ,  _  2S2g3  e2 

0 

1 

-1 

£2  +  £3  h 

(0  -1 

1  / 

The  stability  of  the  different  fixed  points  in  the  £1,2,3 
parameter  space  is  inferred  from  the  scaling  dimensions  of 
Table  I,  and  shown  in  Fig.  5.  Since  each  stable  fixed  point 
implies  a  phase  of  the  junction,  we  shall  also  refer  to  the 
graph  defining  the  regions  of  stability  as  a  phase  diagram. 
It  is  convenient  to  illustrate  this  three-dimensional  phase 
diagram  by  some  cross  sections.  As  all  the  scaling  dimensions 
of  leading  irrelevant  operators  have  a  cyclic  symmetry  on 
wire  indices,  the  stability  regions  show  a  threefold  rotation 
symmetry  around  the  £1  =  g2  =  £3  axis.  Thus,  a  natural  choice 
for  these  cross  sections  is  given  by  planes  normal  to  the 
equal-g  axis:  (£1  +  £2  +  #3)/ 3  =  g,  where  the  parameter  g, 
the  average  of  the  three  Luttinger  parameters,  labels  each 
cross  section.  Excluding  negative  gi ,  these  cross  sections  are 
equilateral  triangles  shown  in  Fig.  5. 

We  first  notice  that  the  decoupled  fixed  point  becomes 
predominant  when  £  <  2/3,  in  which  none  of  other  fixed 
points  are  stable.  One  can  also  show  that  chiral  fixed  points 
appear  when  £  >  1  and  the  Dirichlet  fixed  point  only  appears 
when  £  >  3.  From  Figs.  5(b),  5(d),  and  5(f),  we  observe  that 
the  decoupled,  chiral  and  Dirichlet  fixed  points  are  realized 
around  the  equal-g  axis  forg<l,  l<g<3,  and  g  >  3, 
respectively.  These  results  are  consistent  with  a  junction  of 
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three  identical  TLL  wires  (indicated  as  black  points  at  the 
center  of  triangles),  and  show  how  far  the  Luttinger  parameters 
can  deviate  from  equal-g  axis  before  these  phases  break  down. 
We  find  that,  in  most  regions,  these  phases  are  stable  under  a 
small  perturbation  away  from  the  equal -g  line.  This  is  relevant 
for  the  realization  of  these  phases  experimentally,  as  the  TLL 
wires  attached  to  a  junction  are  likely  nonidentical. 

As  shown  in  Figs.  5(c)  and  5(e),  the  asymmetric  fixed  points 
become  important  around  two  marginal  points  g  1,2,3  =  1,3. 
As  a  small  deviation  of  Luttinger  parameters  from  these  two 
points  easily  realizes  and  switches  between  asymmetric  fixed 
points,  the  resultant  fixed  points  are  thus  highly  sensitive  to 
all  the  three  Luttinger  parameters,  and  not  just  the  averaged 
Luttinger  parameter  g.  Therefore,  precise  control  over  the  TLL 
wires  become  essential  around  these  points.  Note  that  it  may  be 
possible  to  alter  between  different  asymmetric  fixed  points,  and 
form  a  nanoswitch  if  one  can  tune  the  Luttinger  parameters. 
We  mention  in  passing  that  the  Luttinger  parameters  of  wires 
can,  in  principle,  be  modified  by  an  external  gate  capacitively 
coupled  to  the  wire.44 

Finally,  it  is  worthwhile  to  compare  our  findings  with 
those  of  Aristov  and  Wolfle, 58,59  where  two  identical  wires 
are  connected  to  a  wire  with  unequal  Luttinger  parameter, 
by  setting  g\  =  g2  =  g.  In  the  repulsive  and  weak  attractive 
interaction  regime,  g  ~  £3  <  3,  in  which  the  N,  x±,  and  Af¬ 
fixed  points  predominate,  our  results  show  excellent  agreement 
with  their  findings.  In  the  strong  attractive  interaction  regime, 
g  ~  #3  >  3,  the  D  fixed  point  identified  in  the  present  work 
and  in  Ref.  27  (by  nonperturbative  boundary-conformal-field- 
theory  methods),  however,  was  not  found  by  the  perturbative 
renormalization  group  approach  of  Refs.  58  and  59. 

VI.  CONDUCTANCE  RENORMALIZATION  FOR  WIRES 
CONTACTED  TO  FERMI-LIQUID  LEADS 

The  linear  conductances  of  different  fixed  points  were 
calculated  in  Secs.  IIIC  and  V.  Here,  we  discuss  the  effect 
of  attaching  the  wires  to  Fermi-liquid  leads.  Remarkably, 
we  find  that  the  conductance  of  each  fixed  point,  in  the 
presence  of  Fermi-liquid  leads,  renormalizes  to  values  that 
are  independent  of  the  Luttinger  parameters.  This  generalizes 
the  following  interesting  effect  for  the  TLL  quantum  wire  with 
Luttinger  parameter  g:  When  attached  to  leads,  the  measured 
conductance  is  quantized  at  e2/h,  which  is  different  from 
ge2/h.  This  discrepancy  has  been  resolved  by  Maslov  and 
Stone  and  by  Safi  and  Schulz  in  Refs.  38  and  39.  There, 
they  studied  an  inhomogeneous  Luttinger  liquid  and  concluded 
that  the  conductance  of  a  TLL  wire  will  only  depend  on  the 
Luttinger  parameter  at  the  contact.  As  a  Fermi-liquid  (metal) 
contact  can  be  thought  of  as  a  TLL  with  Luttinger  parameter 
g  =  1,  the  measured  conductance  becomes  simply  e2/  h. 

An  alternative  way  to  understand  this  renormalization  of 
conductance  due  to  a  Fermi-liquid  contact  is  to  introduce  a  con¬ 
tact  resistance,  l/Gc  =  (g  —  l)(2g)(h/e2),  at  both  contacts 
in  series  with  the  theoretical  predicted  resistance  l/g(h/e2). 
This  would  give  the  total  conductance  exactly  at  e2 / h.  Here, 
the  contact  conductance  causes  a  voltage  drop  when  matching 
with  the  wire  and  gives  rise  to  a  current/voltage  relation 

U  =  G^Vi  -  V-),  (6.1) 


Fermi  liquid  Luttinger  liquid 


FIG.  6.  (Color  online)  The  voltages  for  the  incoming  and  outgo¬ 
ing  chiral  currents  in  the  Fermi-liquid  and  connected  TLL  wire. 

where  Vi  is  the  potential  of  the  carriers  injected  into  wire  i 
and  Vi  is  the  applied  voltage  at  the  contact  connected  to  the 
wire  i.  We  can  understand  this  result  by  thinking  in  terms 
of  the  right-  and  left-moving  currents  inside  the  wire  and  the 
Fermi-liquid  lead  as  seen  in  Fig.  6.  The  incoming  electrons 
from  the  Fermi-liquid  side  are  at  the  voltage  V )  of  the  reservoir. 
The  outgoing  current  may  be  at  a  different  potential  Vtout  at 
the  contact  point,  even  though  the  electrons  are  expected  to 
relax  to  the  equilibrium  voltage  as  they  propagate  in  the  lead. 
Similarly,  the  right-  and  left-moving  currents  in  the  TLL  wire 
are  respectively  at  voltage  Vi  and  V™ .  The  junction  of  the 
Fermi-liquid  lead  and  the  TLL  has  an  effective  conductance 
gee2 /h  with 

\/ge  =  (1  +  l/g)/2. 

The  current  is  then  related  to  the  difference  in  the  voltages  of 
the  two  incoming  currents  as 

h  =  ge(Vi  -  vr). 

For  the  TLL  wire,  we  can  similarly  write 

h = g(Vi  -  vr). 

Combining  the  above  two  equations  leads  to  Eq.  (6.1). 
Generally,  it  is  useful  to  define  a  contact  conductance  tensor 
Gc  =  8ij(e2/ h)(2gi)/(gi  —  1)  for  a  junction  of  three  wires. 

As  the  measured  conductance  is  based  on  the  applied 
voltage  at  contact,  one  can  define  a  renormalized  conductance 
tensor,  Gij ,  as  7/  =  G ij  Y/  •  To  connect  this  conductance 

with  one  we  found  in  the  previous  subsection,  we  invoke 
current  conservation 


h  =  EG‘/V/  =  EFY- 

j  j 

(6.2) 

Together  with  Eq.  (6.1),  one  can  show  that27 

g  =  (1  +  gg;1)_1g, 

(6.3) 

or  equivalently 

G_1  =  G"1  +  g;\ 

(6.4) 

which  has  a  simple  interpretation  of  resistances  connected  in 
series. 

Besides  the  decoupled  fixed  point  that  has  obvious  vanish¬ 
ing  conductance  GN  =  0,  we  shall  now  apply  Eq.  (6.3)  and 
obtain  the  measured  conductances  of  the  other  fixed  points.  For 
Dirichlet  fixed  point,  we  have  the  renormalized  conductance 

G"  =  jjdS.j  ~  1).  (6.5) 

As  the  largest  conductance  one  can  obtain  from  single¬ 
particle  unitary  scattering  is  GUjk  =  —(4/9 ){e2  /  h)  for  j  /  k,20 
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the  enhanced  conductance  above  demonstrates  the  role  of 
multiparticle  scattering  processes. 

Upon  attaching  the  wire  to  external  Fermi-liquid  leads,  the 
measured  conductances  of  chiral  fixed  points  become 

1  e2 

G%  =  (6.6) 

(i.e.,  the  currents  flow  only  from  lead  1  to  2,  2  to  3,  and  3  to 
1  for  the  x+  fixed  point  and  in  reversed  order  for  x~  fixed 
point).  Finally,  the  measured  conductance  of  asymmetric- At 
fixed  points  reads 

Gfj  =  j(~l  +  Sij  +  Sik  +  2 Sjk  -  3 SijSik),  (6.7) 

which  simply  indicates  a  decoupled  wire  i  with  the  rest  of  two 
wires  fully  conducting. 

We  note  that  all  renormalized  conductance  GD ,  Gx±,  and 
GAi ,  are  the  same  as  the  unrenormalized  conductance  G  with 
all  gi  =  1.  This  result  highlights  that  the  dc  conductance  of 
a  junction  of  TLL  wires  depends  only  on  the  asymptotic 
value  of  the  Luttinger  parameters  of  the  wires,  and  in  the 
case  when  Fermi-liquid  leads  are  attached,  this  asymptotic 
value  is  gt  =  1.  Thus,  when  in  contact  to  Fermi-liquid  leads, 
the  conductance  tensor  of  the  junction  will  take  on  the  universal 
values  listed  above  for  the  different  fixed  points.  Notice, 
however,  that  which  fixed  point  is  selected  still  depends  on 
the  gt  of  the  interacting  wire  segments. 

VII.  CONCLUSION 

In  summary,  we  applied  the  DEBC  method  to  a  two- wire 
and  a  Y  junction  of  TLL  wires  with  generally  unequal 
Luttinger  parameters.  Lor  two  spinless  wires,  we  successfully 
reproduced  the  prediction  that  all  properties  of  the  junction 
are  determined  by  a  single  effective  Luttinger  parameter  ge. 
We  verified  this  prediction  by  direct  numerical  calculations  on 
the  lattice  with  the  method  of  Refs.  66  and  67:  We  observed 
numerically  that  as  long  as  ge  <  1 ,  even  if  one  of  the  wires  has 
attractive  interactions  gi  >  1 ,  any  impurity  leads  to  a  vanishing 
linear  conductance.  Moreover,  we  found  that  the  nonuniversal 
corrections  to  the  correlations  across  the  junction,  which 
come  from  perturbations  with  irrelevant  boundary  operators 
to  the  decoupled  fixed  point,  are  independent  of  the  individual 
Luttinger  parameters  and  only  depend  on  ge  and  the  local 
microscopic  structure  of  the  junction.  Lor  ge  >  1 ,  we  explicitly 
found  a  universal  conductance  of  gee2 /  h  regardless  of  the 
individual  Luttinger  parameters  and  the  microscopic  details. 

Lor  a  Y  junction  of  nonidentical  TTLs,  we  found  that  the 
N-BC,  x_BC,  and  the  D-BC  are  stable  within  regions  of  the 
(gi,g2,g3)  parameter  space,  which  we  explicitly  determined. 
By  identifying  three  more  asymmetric  fixed  points,  corre¬ 
sponding  to  only  one  decoupled  wire,  and  determining  the 
region  of  stability  of  each  fixed  point,  we  determined  the  full 
phase  diagram  of  the  Y  junction.  We  also  obtained  explicit 
formulas  for  the  conductance  of  all  fixed  points.  The  findings  of 
this  work  have  direct  experimental  relevance.  In  particular,  our 
results  shed  light  on  the  issue  of  connecting  interacting  TLLs  to 
Lermi-liquid  leads  for  the  measurement  of  transport  properties. 
Our  work  also  provides  an  important  theoretical  extension  of 


the  well-known  results  on  transport  through  junctions  of  two 
and  three  identical  TLL  wires. 

As  an  outlook,  we  finally  discuss  how  to  generalize  the 
current  method  to  a  junction  of  TV  >  3  wires.  The  key  to  such 
generalization  is  to  identify  the  possible  fixed  points  via  the 
rotation  matrix  1Z.  These  rotation  matrices  are  constrained 
by  charge  conservation  and  can,  in  general,  be  expressed  as 
SO(N  —  1)  matrices  after  eliminating  the  total-charge  (center- 
of-mass)  mode  (cf.  Appendix).  Then,  one  has  to  identify  the 
corresponding  N-BC,  D-BC,  chiral-like  BC,  and  asymmetric 
BC  and  analyze  their  stability.  One  useful  trick  for  identifying 
the  possible  stable  fixed  points  is  to  utilize  the  fact  that  a  stable 
BC  would  make  certain  boundary  operators  effectively  equal 
to  identity  (cf.  the  discussion  in  Appendix  A3).  Of  course, 
with  the  number  of  wires  increasing,  the  number  of  possible 
fixed  points  also  increases  and  the  analysis  becomes  more 
complicated. 
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APPENDIX:  JUNCTION  OF  THREE  WIRES  WITH 

UNEQUAL  LUTTINGER  PARAMETERS:  DELAYED 
EVALUATION  OF  BOUNDARY  CONDITIONS  ANALYSIS 

In  this  Appendix,  we  present  the  details  of  the  DEBC 
analysis  of  the  stability  of  the  following  fixed  points  for 
a  junction  of  three  quantum  wires:  decoupled,  chiral-x±, 
Dirichlet,  and  asymmetric-A;  fixed  points.  The  system  is 
described  by  the  action  (2.1)  with  i  =  1,2,3  and  the  hopping 
Hamiltonian  (2.10)  with  atj  =  y  /3,  where  y  is  the  magnetic 
flux  through  the  ring  at  the  junction. 

To  simplify  the  notation,  we  drop  the  overhead  tilde  (~) 
symbol  for  the  rescaled  fields  throughout  this  Appendix.  (All 
the  bosonic  fields  are  rescaled.)  To  employ  the  DEBC  method, 
it  is  convenient  to  choose  a  proper  basis  for  the  rescaled  boson 
fields.  In  the  first  step,  we  identify  the  following  center-of- 
mass  field,  which  always  satisfies  the  N-BC  due  to  charge 
conservation: 

$o  =  /  ,  - - (VgTffi  +  \fg2V2  +  Vg3<P3)-  (Al) 

V£  i  +  82  +  #3 

The  dual  field  to  d>o,  i.e., 

®o  =  (\fg\Q\  +  \[g2®2  +  \fg3Q3)  (A2) 

V#  1  +  #2  +  #3 

then  becomes  a  constant  and  can  be  simply  neglected.  We  then 
define  another  two  orthonormal  boson  fields, 

.  _  sf&m  -  sfiivi 
1_  VsTTgl  ’ 
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,  ly/gigi<P\  +  ^fglgi^l  -  (g  1  +  gl)(P3]  ,  A  ry x 

^2  =  - ,  - ,  (A3) 

V#1  +  82  +  g3\/gl  +  g2 

as  well  as  their  dual  fields,  ©i;2-  Note  that  the  choice  of 
basis  above  is  arbitrary  but,  as  will  become  apparent,  is  a 
convenient  one.  We  will  organize  the  fields  above  into  a  vector 
4>  =  (Oi,(b2)r  and  its  dual  vector  0  =  (©i,©2)r. 

It  is  useful  to  define  the  following  vectors: 


K  i 

K2 

k3 


VgT  Vgl  +  g2  +  g3  ^ 
\/2g2(gl  +  g2)  ’  \/2g3(gl  +g2)/’ 

Vg2  Vgl  +  g2  +  g3^ 
V2gl(gl  +g2)’  V2g3(gl  +g2>/’ 


Vgl  +  g2 
V2gigi 


,0 


(A4) 


which  will  further  simplify  the  notation.  Notice  that  these  three 
Kj  vectors  add  up  to  0  for  any  gt . 

By  using  the  notation  in  Eq.  (3.3)  and  neglecting  ©o  (due 
to  the  N-BC  on  O0),  the  boundary  operators  for  single-particle 
processes  are  categorized  in  four  classes,  and  listed  in  Table  II. 
As  in  the  two-wire  case,  the  higher-order  processes  can  be 
constructed  from  these  single-particle  boundary  operators. 

Using  the  ansatz  (3.2),  we  then  write  a  rotation  matrix 


n,= 


(cos§ 
-sin£ 


(A5) 


where  §  is  a  rotation  angle.  The  rotation  matrix  above  relates 
0L  =  (O  +  0)/2  to  (j)R  =  (4>  —  ©)/2.  In  terms  of  7 Z%,  the 
scaling  dimension  of  boundary  operators  is  given  by  Eq.  (3.5). 
We  now  proceed  to  the  stability  analysis  of  the  four  Y-junction 
fixed  points. 


1.  Decoupled  fixed  point 

The  decoupled  fixed  point  corresponds  to  the  N-BC  for  <I> 
field,  and  makes  the  0  field  a  pure  number.  Therefore,  all 
backscattering  processes  are  effectively  identity.  The  rotation 


matrix  is  simply  equal  to  7 ^=0  =  1.  From  Eq.  (3.5),  the  scaling 
dimension  of  an  arbitrary  operator  with  the  N-BC  becomes 
A@b  =  \a\ 2 .  The  explicit  scaling  dimensions  for  the  operators 
in  Table  II  are  listed  below. 


Ob 

Ags  (N-BC) 

tRL  tRL  tLL  tRR 

121  ’M2  ’  Ml  ’  Ml 

IM2  =  (gi+g2)/(2g1g2) 

r rRL  tRL  tLL  tRR 

M2  ’  M3  ’  M2  ’  M2 

1^1  I2  —  (g2  +  gl)/Og2gl) 

r rRL  ’-rRL  tLL  tRR 

M3  ’  Ml  ’  1 13  ’  1 13 

IM2  =  (g3  +  gi)/(2g3gi) 

Note  that  condition  Aqb  >  \  determines  the  stability  region 
of  the  N-BC,  shown  in  red  in  Fig.  5. 


2.  Dirichlet  fixed  point 

The  Dirichlet  fixed  point  corresponds  the  D-BC  on  the  O 
field  (i.e.,  4>  is  effectively  a  constant  at  boundary).  The  rotation 
matrix  of  D-BC  simply  reads  7 Zj?=n  =  —1.  From  Eq.  (3.5),  the 
scaling  dimension  of  an  arbitrary  operator  with  D-BC  becomes 
Aqb  =  \b\2.  Unlike  the  N-BC,  none  of  the  single-particle 
processes  becomes  identity  with  D-BC.  However,  some  of 
the  two-or-more-particle  processes  do  become  identity  under 
D-BC,  for  instance:  TRLTR2L\  T32T23\ and  TrlTrl] .  This 
indicates  that  the  Dirichlet  fixed  point  is  associated  with  the 
Andreev  reflection. 

Hereby,  we  only  list  the  scaling  dimensions  for  the  =b  cycle 
as  they  are  the  leading  irrelevant  operators. 


Ob 

(D-BC) 

tRL  tRL 

12\  ’M2 

g3(gl  +  gl)/2(gl  +  g2  +  g2) 

' -rRL  'T'RL 

M2  1  1 22 

81(82  +  83) /^(8i  +  82  +  £3) 

T'RL  tRL 

1 13  ’  Ml 

g2(g3  +  gl)/2(gl  +  g2  +  g3) 

TABLE  II.  The  boundary  operators  corresponding  to  single¬ 
particle  processes  at  the  Y  junction. 


(a)  =b  cycle 


VST8283  /s...  ^ 

tRL  n 
1  21(12) 

-  e±iK*-*e 

VS 18283  via 

tRL  n 
1 32(23) 

-  e±iK^e 

g 81+82+83  ^ 

VS 18283  /a 

tRL 

1 13(31) 

.  e±iK 2*J 

Vsl+82+83  {Z  2) 

(b)  1 


tRL  n 
1 11 

—i 

~  e 

:  zV8]8283 

'  a/ 81+82+83 

(zxK  i>0 

:  2V8]8283 

(zxK2)-& 

tRL  r- 
1  22 

~  e 

J 81+82+83 

,  2V8\8283 

(zxK3)-0 

tRL  n 
1  33 

~  e 

y/8\+82+S3 

(c)  LL-RR  processes 

f/pm  '  +j(yAsi^&l+  Vg3(gl-g2)®2  s 

T^L(KK)  ^  giKi  Vsi+82  y/2(gi  +g2)(gl  +g2+g3) ’ 

_  ygmQi  yg3(gi+2g2)Q2  a 

Tyf{RR)  ~  eiKv*e  W%g\+Sl)  V2(gl+g2)(gl+g2+g3) 

_  Vglg20l  ,  Vg3(2gl+g2)®2  ^ 

TLL{rR)  ^  eiK2  ^e+  V2(gl+g2)  V2(gl+g2Xgl+g2+g3) 


Now,  the  Dirichlet  fixed  point  is  stable  only  when  all  these 
scaling  dimensions  A®b  >  1.  In  Fig.  5,  the  stability  region  of 
Dirichlet  fixed  point  (D-BC)  is  painted  in  green. 


3.  Chiral-x±  fixed  points 

The  chiral- x±  fixed  points  are  defined  as  follows:  the  BC 
corresponding  to  the  x  +  (x)  fixed  point  would  effectively 
make  all  operators  in  +  (— )  cycle  equal  to  identity  for  all 
Luttinger  parameters.  With  this  in  mind,  we  could  derive 
the  following  relationship  for  the  rotation  angles  of  the 
corresponding  rotation  matrices: 


tan§±  =  d= 


/  gl  +  g2  +  #3 


glg2g3 


(A6) 


where  £±  are  the  rotation  angles  of  x±-BC,  respectively.  The 
rotation  matrices  are  obtained  by  plugging  in  the  respective 
rotation  angles  into  Eq.  (A5).  Here,  we  do  not  show  them 
explicitly. 
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Let  us  first  focus  on  the  x+ 

fixed  points.  The  scaling 

dimensions  of  the  single-particle 
cycle)  are  as  follows. 

processes  (excluding  + 

Ob 

TRL  rrRL  'T'LL  ’-rRR 

Ml  ’  123  »  1  13  ’  Ml 

Zgl(g2+g3) 

glg2g3+gl+g2+g3 

TRL  tRL  7 -’LL  ’t’RR 

2g2(gl+g3) 

1 22  »  *31  »  Ml  ’M2 

glg283  +8l+82+83 

’T’RL  ’T’RL  j^LL  t’RR 

Zg3(gl+g2 ) 

1  33  ’  1 12  ’  1  31  ’  1  13 

g\g2g3+g\+g2+g3 

Notice  that  these  scaling  dimensions  are  cyclic  in  three  indices 
and  hence  all  operators  are  important  for  determining  the 
stability  of  the  x+  fixed  point. 

As  for  the  X-  fixed  point,  one  can  show  that  all  the  leading- 
order  operators  have  exactly  the  same  scaling  dimensions  listed 
in  the  table  above.  Thus,  both  x±  fixed  points  share  exactly 
the  same  stability.  In  Fig.  5,  the  stability  region  of  chiral- x± 
fixed  point  is  shown  in  orange. 

4.  Asymmetric  fixed  points 

Although  the  use  of  rotation  matrices  is  useful  for  other 
fixed  points,  it  is  most  convenient  to  identify  rotation  matrices 
directly  in  the  rescaled  boson  field,  4>R,R,  basis.  Because 
the  decoupled  wire  effectively  has  the  N-BC  for  itself  and 
the  connected  wires  should  follow  mutual  D-BC,  by  using 
Eq.  (3.14),  the  rotation  matrix  of  asymmetric  fixed  point,  A\, 


has  the  form 

n  °  0  \ 

TLa'  = 

Q  82~  g3  2  f 8283 

82+83  82+83 

Q  V#2#3  g3-g2 

\  82+83  82+83  / 

(A7) 

while  those  of  A2,3  can  be  constructed  by  permuting  the  indices 
in  the  corresponding  matrix  elements. 

By  using  this  rotation  matrix  (A7),  it  is  straightforward  to 
show  that  the  following  single-particle  tunneling  processes  are 
equal  to  identity:  T^L,  T^L,  and  T^L.  In  addition,  the  scaling 
dimensions  of  the  leading  relevant/irrelevant  operators  read  as 
follows. 

Ob 

A A0\  (ArBC) 

’T’RL  ’T’RL  T’LL  t’RR  'T'RL  'T'RL  ’T’LL  t'RR 

1 21  ’  1 12  ’Ml  ’Ml  ’  1 13  ’Ml  ’  1 13  »  1 13 

r rRL  ’-rRL  rrLL  'yRR 

M 2  ’M3  ’M2  ’  *32 

rrRLrrRLJ\  rrRLrrRLJ\ 
l2\  1  12  ’  1 13  J31 

gl+g2+g3+glg2g3 

Zgl(g2+g3 ) 
2g283 
g2+g3 

2(3-  +  — f — ) 

vg  1  g2+g3  7 8 9 10 * 12 13 

Here,  we  notice  that  some  leading  order  operators  are 
two  particle  processes.  To  obtain  the  scaling  behaviors 
of  operators  near  the  A2j 3  fixed  points,  one  can  simply 
permute  the  indices  of  the  Luttinger  parameters  with  the 
corresponding  operators.  In  Fig.  5,  the  stability  regions  of 
A  1,2,3  fixed  points  are  shown  in  yellow,  gray,  and  blue, 
respectively. 
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